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Abstract: Based on the discussion of relationships between thinning and wind damage, and published information, a method for 
estimating risk ratios of wind damage was developed. Estimations of risk-ratio for Pinus thunbergii trees and stands were de¬ 
duced from stem bending theory and coefficients characterizing wind profile, distribution of branches and optical stratification po¬ 
rosity. The results showed that if the value of constant f) in the branch distribution-model equals the attenuation coefficient a 5 in 
the wind profile model for a single tree crown, then the parameter H/Di 3 3 (height over stem diameter cubed) can be used to 
compare and evaluate the risk-ratio of wind damage for individual trees. The same method can be applied to stands using the 
coefficient of wind profile in a stand, i.e. attenuation coefficient a , the coefficient from distributions of optical stratification porosity, 
i.e. extinction coefficient v , and the parameter Di 3 3 . The application of parameter H/Di 3 3 and the process of determining risk ra¬ 
tios of wind damage for stands were also given in the paper. 
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Introduction 

Wind damage is one of the major constraints on the 
practice of silviculture (Blackburn and Petty 1988a; Mayer 
1989; Rollinson 1989; Quine 1995; Pellola et al. 2000). 
How wind interacts with forests concerns the forest man¬ 
ager in a number of important ways. The influence of wind 
on forest management and planning can be considered at 
two main stages in the life of the forest: during and after 
regeneration (Godwin 1968). Obviously, the influence of 
wind on the latter, i.e., established stands, has been paid 
great emphasis, and many studies have been reported 
(Hutte 1968; Stumbles 1968; Petty and Worrell 1981; 
Cremer et al. 1982; Petty and Swain 1985; Galinski 1989; 
Peltola and Kellamaki 1993; Quine 1995; Gardiner et al. 
1997; Peltola et al. 1999; Moore and Quine 2000). Almost 
all of the studies are concerned with timber-production 
forests. Forest managers, however, need to understand 
effects of wind in forests used for non-timber purposes. 
Coastal forests, for example, alter the wind thereby pro¬ 
viding protection as well as many other benefits such as 
recreation, nature conservation, timber, fungi, and berries. 
Therefore, it is desirable and necessary to make coastal 
forests resistant to severe damages, especially those re- 
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suiting from strong wind. However, research in this regard 
is poorly developed because it has been thought that there 
is no countermeasure against meteorological extremes 
that are beyond human control (Matsuzaki 1994). Thinning 
is one of the few countermeasures against wind damage 
that can be considered in the management of a forest. 
However, outcomes of thinning in coastal forest are difficult 
to determine because, on one hand, thinning can improve 
the stability of the plantation forest, whiie on the other hand 
risk of wind damage increases immediately after any thin¬ 
ning (Cremer et al. 1982; Gardiner et al. 1997), especially 
in the coastal area where forests suffer from exposure to 
relentless winds blowing off the sea (Perry 1994; Zhu et al. 
2000b). 

Although thinning may result in high probability of wind 
damage, it is necessary to have narrow initial spacing for 
the survival and growth of the trees in the ruthless coastal 
area. Therefore, it is important to understand the relation¬ 
ships between thinning and wind damage in detail, and 
recognize how to estimate the risk of wind damage for a 
coastal forest stand with various thinning intensities, and 
what damage will be caused by extreme winds. Managers 
need basic information on the influences of wind to solve 
the problems of establishment and management for the 
coastal forests. Such information could also be useful for 
researchers of micro and meso-scale meteorology, town 
and regional planners, and farmers. An understanding of 
the effects of wind on forests may also be useful for as¬ 
sessing forest health and stability. 

We recently completed a study of the influence of thin¬ 
ning on the incidence of wind damage in a coastal pine 
(Pinus thunbergii Pari.) plantation. The objective of the 
study is to identify (1) how wind damage is related to thin- 
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ning, (2) what factors influence the stability of individual 
trees and stands, (3) what degrees of risk might be quanti¬ 
fied after thinning, and (4) how thinning strategy should be 
designed to reduce risks of wind damage for the manage¬ 
ment of coastal pine forests. Here we present a review of 
relationships between wind damage and stand density or 
thinning in timber-production forests, and develop a 
method of wind damage-risk estimation for individual trees 
and stands on the basis of published information (Black¬ 
burn and Petty 1988b; Galinski 1989). 

Thinning and the risk of wind damage 

Wind damage can be mainly classified into stem dam¬ 
age and canopy damage. The most common report is stem 
damage, which includes uprooting, breakage, and bending 
or leaning (Everham 1995). Canopy damage may be quan¬ 
tified directly on the basis of branch loss, canopy damage 
or defoliation. 

Generally, thinning entails risk of wind damage, i.e., the 
likelihood of damage due to strong winds, in the short term 
(Cremer et at. 1982; Quine 1995). Particularly, it could 
have an important and unanticipated impact on the occur¬ 
rence of wind damage where the climate is windy and sites 
are poor (Quine 1995). It has been reported that the stabili¬ 
ties of trees and stands are reduced after thinning, and the 
degree of the destabilization, i.e., the risk of wind damage 
caused by thinning, depends not only on meteorological 
factors but also on site conditions (topography, soil type, 
soil water content etc.), forest structure (forest age, height, 
density etc.) and the pattern and intensity of thinning 
(Cremer et al. 1982). The relationships between thinning 
and the risk of wind damage can be summarized from lev¬ 
els of both single tree and stand as follows. 

1) The removal of neighboring trees results in an open¬ 
ing of the canopy, and the wind gusts penetrate deeper. 

2) The wind protection from the neighboring trees and 
the mechanical crown contact are reduced by thinning. The 
two kinds of reductions contribute to increased tree sway 
amplitudes, i.e., the increase of destructive storm moment 
(Gardiner 1995; Nielsen 1995). 

3) It is the thinning that increases the aerodynamic 
roughness of the canopy surface and turbulence. This will 
enhance the momentum transfer from wind to stand. The 
increased amount of transferred energy is distributed to a 
reduced number of trees per hectare, i.e., the overall wind 
load has to be shared by the reduced number of remaining 
trees. Thus, thinning increases the destructive storm mo¬ 
ment of the retained trees (Cremer et al. 1982; Nielsen 
1995; Peltola 1996a). 

It has been suggested that the most dangerous situation 
arises when the wind spectrum includes strong gusts 
whose frequency coincides with the natural frequency of 
oscillation of the trees (Peltola et al. 1993; Gardiner 1995). 
It may indeed be totally inadvisable for thinning in stands 
which are already weak because of overstocking, and 


perhaps also in stands where root development is ex¬ 
tremely restricted (Cremer et al. 1982; Ruel et al. 1998). 

Conversely, thinning reduces the risk of wind damage by 
improving the stand stability in the long term. With the 
growing of trees, the thinned stand will gradually regain the 
stability and become stronger than before the thinning 
(Matsuzaki 1994). The increased diameter growth 
achieved at lower stocking or density is generally higher 
than height growth (Deans and Milne 1999). This attributes 
not only to the fact that each tree gains better access to 
light (Cremer et al. 1982), water and nutrients (Kuraji et al. 
1997), but also to the fact that the increased exposure to 
wind results in more swaying (Cremer et al. 1982; 
Telewsiki 1995). A study of thigmo-morphogenesis 
(Telewsiki and Jaffe 1981) indicated that swaying in¬ 
creased diameter growth and decreased height growth. 
Thus, thinning can reduce the ratio of H/Di 3 or increase 
taper (Di 3 /H) of trees (H is tree height, Di 3 is diameter at 
breast height, 1.3 m). Exposure to wind also tends to in¬ 
crease growth of roots (Grace 1988; Stokes et al. 1995), 
and thus improves the tree strength and stability. Long 
after thinning, stands redevelop a smooth canopy surface, 
low momentum transfer to the canopy, high wind protection 
between the neighboring trees, and high reduction of tree 
sway amplitudes and increment of strength. 

The time needed to regain certain stability for thinned 
stands depends on the growth rate, stand age and the 
state of the stand at the time of thinning. Generally, the 
changes of thinning are the most significant when the dura¬ 
tion of their effect (time for crowns to grow back) is not 
longer than recurrence interval of the appropriate threshold 
wind speed (Quine et al. 1999). 

In a word, thinning is one of the most important coun¬ 
termeasures for the sustainable forestry to prevent wind 
damage (Matsuzaki and Nakata 1993); however, the forest 
soon after thinning is vulnerable to strong wind because it 
does not quickly recover resistance to extreme wind. Re¬ 
solving the dilemma between the benefits and risks of thin¬ 
ning may require management skill, and luck. 

Index for the risk of wind damage 

The risk of wind damage is subject to very complex fac¬ 
tors associated with the stature of trees, the aerodynamic 
properties of a stand, the general and local incidence of 
winds, soil conditions, and so on. These factors can be 
classified as three main aspects, i.e., weather, site and 
stand conditions (Cremer et al. 1982). The weather condi¬ 
tions are always beyond human control (Matsuzaki 1994), 
the site conditions are mainly permanent for a special loca¬ 
tion, while the stand conditions are the most important and 
controllable, and the features of a stand are subject to 
change with growth and silvicultural measures. Therefore, 
the important effects of stand structure or configuration on 
wind loading have been dealt within many studies (Cremer 
et al. 1982; Petty and Swain 1985; Galinski 1989; Gardiner 
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etal. 1997). 

Foresters have used the ratio of H/D1.3 or taper Di 3/H as 
a measure of risks of wind and snow damage (Valinger and 
Pettersson 1996). The ratios of H/D1.3 were calculated in 
several ways, but mostly from the mean height and di¬ 
ameter of all the trees in the stand. However, Cremer et at. 
(1982) suggest that both tree height (H) and diameter (Di 3) 
calculated from the mean values of the 200 largest stems 
per hectare, i.e., H 200 L/D 200 L (where H 2 oot.and D 200 L are the 
mean tree height and mean Di 3 of the largest 200 stems 
per hectare) can be considered as a valuable index of the 
risk of wind damage. This is because that smaller trees are 
far less significant than the dominant trees in determining 
the stability of the stand, and selective removal of the 
smaller trees by thinning will at once boost the ratio of 
Di 3/H, and thus may wrongly suggest that the risk is low¬ 
ered. Ratio of H/D1.3 is considered as the most important 
factor likely to influence the stability of stand (Petty and 
Worrell 1981; Cremer et at. 1982; Blackburn and Petty 
1988b). 

Other ratios of height and diameter have been used as 
alternatives to H/D 13 , for example, H/Di 3 3 , H/D 1 . 3 2 and 
H 3/2 / Di 3. Generally, H/Di 3 3 relates to the strength of a 
cantilever of uniform resistance (Gardiner et at. 1997). 
H/Di 3 Z is closely related to the shape of tree stem (Cremer 
et at. 1982). And H /Di 3 gives the greatest relative weight 
to tree height and relates best to elastic strength properties 
(Valinger etal. 1993). 

The relationships between the index values and inci¬ 
dence of wind damage showed that the combination of 
height and diameter in H/Di 3 was valuable for providing 
the first indication. However, risk of wind damage depends 
not only on the sturdiness of the tree population but also 
on the aerodynamic properties of the stand, as well as on 
wind conditions, therefore, the relationships between the 
above mentioned index values and incidence of wind 
damage can only provide a broad indication of risk. 

Estimation of risk of wind damage 

Foresters need a practical measure to assess the risk of 
wind damage for a particular tree or stand. Blackburn and 
Petty (1988b) and Galinski (1989) suggest a practical 
method for evaluating the risk of wind damage from the 
combination of stem bending theory and wind regime. The 
fundamentals of this method are accepted in this paper 
also, but some details are developed for estimating the risk 
of wind damage for a single tree (Pinus thunbergii Pari.) 
and for a coastal stand of the same tree species. 

Estimation of risk-ratio for individual tree 

The diagnosis of wind risk for individual trees and stands 
requires evaluation of how well trees or stands are accli¬ 
mated to wind load (Mitchell 2000). Firstly, it is assumed 
that a single tree of P. thunbergii can be modeled as an 
elastic beam with constant modulus of elasticity. Wind act¬ 


ing on the tree crown creates a drag force, which when 
acting along the stem produces increasing turning mo¬ 
ments towards the base of the tree (Galinski 1989). These 
moments are resisted by the strength of stem and the 
root-soil system, and the drag force that causes the tree to 
be bent results from wind action (Galinski 1989; Mitchell 
2000). If it is assumed that the main cause of the drag 
force is the presence of needled tree branches which can 
be approximately characterized by their length (Galinski 
1989), then the momentum exchange between moving air 
(wind) and the tree branches causes the drag force which 
for a given tree is nearly proportional to the wind speed 
(Peltola 1996b). Therefore, the drag force in a single whorl 
or verticil of branches in a single tree can be represented 
by equation (1) (Galinski 1989). 

F,=kLiU z ( 1 ) 

where k is constant coefficient, F, is the drag force (N), L, is 
the total length of needled branches in ilh whorl from the 
base, Uz is wind speed within a single crown correspond¬ 
ing to the height of the /th whorl (m - s' 1 ). 

The length of needled branches is summed over the 
whorl because the influence of the wind direction is not 
considered. Galinski (1989) suggests the principle of equa¬ 
tion (1) for conifer trees is as follows: the whorl is the 
smallest reasonable unit for which the drag force can be 
calculated without considering the behavior of each branch 
under wind action. The leeward branches are bent into a 
parallel position to the wind speed vector, and the wind¬ 
ward branches are strongly deflected from the position by 
wind. Thus, it is reasonable to assume that the whole whorl 
resistance results from wind speed and the sum of needled 
branches in the whorl. 

Wind speed within a single tree crown of P thunbergii 
obeys the exponential form given in following equation 
(Zhu et al. 2000a). 

U 5z =U out exp[-a s (\-H i /H)] (2) 

where H, is the height of interest within the canopy (m), H 
is maximum height of the tree crown (m), U sz is wind speed 
(m-s' 1 ) inside the canopy at height z, U 0 ut is wind speed 
(m ■ s' 1 ) outside the canopy at a certain height (generally at 
2 m, i.e. U ou p U z =2m), it can also be replaced by Uh (wind 
speed at the tree top, m-s' 1 ) because UJU ■ = (1/ 
k )ln[(z-d)/Zo] (logarithm law, U■ is friction velocity, m-s" 1 , 
k is the von Karman constant, d is zero-plane displace¬ 
ment, m, Zo is roughness parameter, m), and a s js a con¬ 
stant, which determines the form of wind profile within the 
single tree crown. 

Replacing U sz in equation (1) with equation (2) results 
in: 

F i = kLjU out expha s (1 - //, / H)] ( 3 ) 

The bending moment in basal cross-section (M,) result¬ 
ing from action of the drag force F, equals: 



312 


ZHU Jiao-jun et al. 


M i = H i F i = kH i L i U oui ex P[-«s(1 - H i / H )] (4) 

and the total bending moment ( M ts ) equals the sum over all 
whorls along the whole tree height. 

=Y, H ‘ Fi= W «“' S H • L ‘ exp[_ “ s (1 " H ‘' H)] ( 5 ) 

However, the parameter of L, is unique for each tree 
species, for P. thunbergii , the distribution of L, increases 
exponentially with decreasing height (Zhu et al. 2000a), 
i.e., 

L, =exp[/3(l - Hj / H)] (6) 

where /3 is constant coefficient. 

Replacing L, in equation (5) with equation (6), equation 
(5) becomes: 

M„ = W oul ^ H t exp[0(l - H, / H)] expl-a, (1 - H, / H)] = 
k U out £ Hi exp [(P - cc s )(1 ■- Hi /H)] (7) 

According to the proposition by Galinski (1989), if every 
transverse cross-section of the stem is circular, the strain 
in its outer wood is given by equation (8) (Galinski 1989; 
Blackburn and Petty 1988a, b). 

S«) = 4M ft /D? 3 * (8) 

where S(t) is strain in stem outer wood (N), ir is Pythago¬ 
ras constant, 3.1416. 

Integrating equations (7) and (8) results in: 

SO) = {4k U out ^ H, exp[(/3 - a s )(l - H, / H)]) / D? 3 * (9) 

For a given wind speed Uout, the differences in tree ar¬ 
chitecture cause a difference in risk of wind damage for 


particular trees. Thus, it is possible to define the risk of 
wind damage as the individual tree architecture contribu¬ 
tion to the strain from equation (9). 

As 4klW * is constant for a given wind speed outside 
the tree crown, the risk-ratio R(s) suggested by Galinski 
(1989) can be defined as, 

R{s) = H, exp[(/3 - a s )(t - H, / H)]} / D ? 3 (10) 

Galinski (1989) suggests that the variation of individual 
ratio, R(s) is relatively low for most trees. This suggests 
that a unique value of R(s) exists for each tree and it is 
independent of the site. Therefore, equation (10) can be 
used to evaluate or compare the risk of wind damage for 
individual trees. 

For tree species of P. thunbergii, we found that if U ou t in 
the model of wind speed within the single-tree crown 
(equation 2) is replaced by Uh, then the value of coefficient 
a s in equation (2) nearly equals coefficient P in equation 
(6). In order to simplify the expression of the problem 
(equation 10), we can assume that coefficient j3 equals 
a s in value, which means that equation (10) becomes: 

R(s) = H/ nli ( 11 ) 

The simple expression of risk estimation of R(s) (equa¬ 
tion 11) allows a comparison of growth strategies among 
individual pine trees. The higher the R(s) is, the higher the 
probability that the considered tree will be destroyed during 
a strong wind. 

The risk estimation of R(s) for the P. thunbergii trees in a 
coastal forest, which was thinned in four treatments, i.e., 0, 
20%, 30% and 50% thinned, which are referred to treat¬ 
ment 1, treatment 2, treatment 3, and treatment 4, respec¬ 
tively, was conducted. Figures 1A and IB are the distribu¬ 
tion of stems by risk ratios of individual trees soon after the 
thinning (December of 1997) and two years after thinning 
(January of 2000). 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
Risk ratio of wind damage 
■® — T reatment 1 —•— T reatment 2 

"*— Treatment 3 * Treatment 4 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 
Risk ratio of wind damage 


Treatment 1 
Treatment 3 


Treatment 2 
Treatment 4 


Fig. 1 The distribution of stems sorted according to risk ratios 

A: soon after the thinning, and B: two years after thinning 
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The proportion of stems with higher risk-ratios of wind 
damage (H/Di 3 3 more than 2.5) plotted in Figures 2A and 
2B showed that stems of higher risk-ratios increased in the 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Risk ratios 

■ Dec of 1997 □ Jan of 2000 


unthinned treatment, but decreased in the 20% thinned 
treatment (the same as in other thinned treatments). 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Risk ratios 

■ Dec of 1997 □ Jan of 2000 


Fig. 2 Changes of the risk ratios 

A: unthinned, and B: 20% thinned 


Estimation of risk-ratio for stand 

Applying the same approach used for a single tree, we 
can assess the risk of wind damage for a forest stand as 
follows. 

If the dynamic wind load is the presence of area of the 
canopy and stem, the momentum exchange between 
moving air (wind) and the tree elements in the stand, which 
causes the drag force, can be estimated from the wind 
speed profiles within the canopy (Peltola et at. 1993; Pel- 
tola 1996b), i.e., 

F d (z) = kF s (z) = k(0.5C d p A z u} z ) (12) 

where Fd(z) is the drag force or dynamic wind load (N), 
F s (z) is the static wind load (N), Cd is the drag coefficient 
(dimensionless), here assumed to be constant, p is the 
air density (kg-m‘ 3 ), A z is the projected area of the crown 
and stem (m 2 ) at height z, and can be approximated as 
equation (13), Ut z is wind speed (m-s' 1 ) at height z within 
the canopy. 

A z =B.dz (13) 

where B z represents the leaf area index (dimensionless). 
As parameters k, Cd and p are constant, let 

k(0.5C d p )= £ (14) 

Then equation (12) becomes, 

F d (z)=SUT z B z dz (15) 

Using the same principle as for the single-tree, the total 
bending moment for a stand, Mu equals: 


H 

zF d (z)dz = £ z £-U } z dz (16) 

) Jo 

Wind speed within the canopy obeys the exponential 
form given in equation (17) (Bergen 1971; Landsberg and 
James 1971; Amiro 1990; Wenzal et at. 1997; Zhu et at. 
2001 a, b). 

U tz =U H exp[-o(l-z/H)] (17) 

where z is the interest height within the canopy (m), H is 
height of top of tree canopy (m), Uh is wind speed (m-s' 1 ) 
at height H, coefficient a is called as attenuation coeffi¬ 
cient. 

B z can be represented by the distribution of optical strati¬ 
fication porosity (OSP) (Zhu et at. 2000b, 2002). 

B z =d{expf-v(l-z/H)]/dz=(v/H)exp[-v(l-z/H)] (18) 

where v is extinction coefficient of distribution of OSP 
(Zhu et at. 2000b, 2002). 

Combining equations (16) and (18) results in: 

J .h 

z exp[-v(l - z /H)]{exp[-o(l - z / H)] } 2 dz 
«0 

f" 

= (eu&v/H) z exp[-(v + 2a)(l - z /H)]dc (19) 

Jtf 0 

where Ho is bole height (m). 

For a given stand, the strain, S(ts) in equation (8) for 
stands, is given by the same form as that of the single tree. 
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S(?j) = (£Uh v/ H) z exp[-(v + 2a)(l- zl H)]ct/ Dj* 3 k (20) 
Jh 0 

As £ Uh 2 v /H jc is constant for a given wind speed at 
the canopy top, the risk-ratio R(t) for a stand can also be 
expressed as 

e» 

R(t)= z exp[-(v + 2 a)(l - z / H)]dz / Dj* 3 (21) 

Jh 0 

D 13 in equation (21) is the mean diameter of the stand 
calculated from various sample patterns (mean Di 3 of all 
stems, mean Di 3 calculated from the mean values of the 
200 largest stems per hectare and so on). Risk estimation 
of R(t) allows a comparison of growth strategies among 
different stands. The higher the R(t) is, the higher the 
probability that the considered stand will be destroyed by 
the strong wind, i.e., R(t) can be used to measure and 
compare the risk of the growth strategy and the thinning 
measures adopted by the stand. The following scheme 
(Fig. 3) presented the process of determination for R(t). 


distribution of optical stratification porosity, combined with 
Di 3 3 could be used to compare the risk of wind damage 
among forest stands. The basic application of the methods 
is to determine the relative risk of wind damage for indi¬ 
vidual trees and stands; in addition, it can be used to 
compare the risk of the growth strategy and the thinning 
measures adopted by the stand. Although the results dis¬ 
cussed above may not be directly applicable to other type 
of forests and other tree species, the risk-ratios of wind 
damage can also be obtained according to special tree 
species using the procedures discussed in this paper. 
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